library(haven)

# Load in Data set
impdat20_relig <- read_spss("./data//Religion_IAT.public.2020.sav")
# Subset to USA
impdat20_relig <- subset(impdat20_relig, countryres_num == 1)

# Meta Info ---------------------------------------------------------------
impdat20_relig$date_mdy <- format(impdat20_relig$date, "%Y-%m-%d")
impdat20_relig$date_mdy <- as.Date(impdat20_relig$date_mdy)
impdat20_relig$week <- format(as.Date(impdat20_relig$date), "%W")

# Demographics ---------------------------------------------------------------
# Age
impdat20_relig$age[which(impdat20_relig$birthyear > 1900)] <- impdat20_relig$year[which(impdat20_relig$birthyear > 1900)] - impdat20_relig$birthyear[which(impdat20_relig$birthyear > 1900)]
impdat20_relig$age_sc <- (impdat20_relig$age-min(impdat20_relig$age, na.rm = T))/(max(impdat20_relig$age, na.rm = T)-min(impdat20_relig$age, na.rm = T))

impdat20_relig$age_cat4 <- NA
impdat20_relig$age_cat4[which(impdat20_relig$age >= 18 & impdat20_relig$age <= 29)] <- "18-29"
impdat20_relig$age_cat4[which(impdat20_relig$age >= 30 & impdat20_relig$age <= 44)] <- "30-44"
impdat20_relig$age_cat4[which(impdat20_relig$age >= 45 & impdat20_relig$age <= 64)] <- "45-64"
impdat20_relig$age_cat4[which(impdat20_relig$age >= 65)] <- "65+"
impdat20_relig$age_cat4 <- as.factor(impdat20_relig$age_cat4)

# Sex
impdat20_relig$sex <- NA
impdat20_relig$sex[which(impdat20_relig$birthSex == 2)] <- "f"
impdat20_relig$sex[which(impdat20_relig$birthSex == 1)] <- "m"

# Race and Ethnicity
impdat20_relig$race <- factor(impdat20_relig$raceomb002, 
                           levels = c(1:max(impdat20_relig$raceomb002, na.rm = T)),
                           labels = c("American Indian/Alaskan Native",
                                      "East Asian",
                                      "South Asian",
                                      "Native Hawaiian/Pacific Islander",
                                      "Black",
                                      "White",
                                      "Other", 
                                      "Multiracial"))

impdat20_relig$hisp <- NA
impdat20_relig$hisp[which(impdat20_relig$ethnicityomb == 1)] <- 1
impdat20_relig$hisp[which(impdat20_relig$ethnicityomb == 2)] <- 0 # drops unknowns

impdat20_relig$race5 <- NA
impdat20_relig$race5[which(impdat20_relig$race %in% c("East Asian", "South Asian", "Native Hawaiian/Pacific Islander") & impdat20_relig$hisp == 0)] <- "Asian/Pacific Islander"
impdat20_relig$race5[which(impdat20_relig$race == "Black" & impdat20_relig$hisp == 0)] <- "Black"
impdat20_relig$race5[which(impdat20_relig$race == "White" & impdat20_relig$hisp == 0)] <- "White"
impdat20_relig$race5[which(impdat20_relig$race %in% c("Other", "Multiracial") & impdat20_relig$hisp == 0)] <- "Other"
impdat20_relig$race5[which(impdat20_relig$hisp == 1)] <- "Hispanic"
impdat20_relig$race5 <- as.factor(impdat20_relig$race5)
impdat20_relig$race5 <- factor(impdat20_relig$race5, levels(impdat20_relig$race5)[c(1:3,5,4)])

impdat20_relig$race3 <- NA
impdat20_relig$race3[which(impdat20_relig$race5 == "White")] <- "White"
impdat20_relig$race3[which(impdat20_relig$race5 == "Black")] <- "Black"
impdat20_relig$race3[which(impdat20_relig$race5 != "Black" & impdat20_relig$race5 != "White")] <- "Other"

# Ideology
impdat20_relig$ideo3 <- NA
impdat20_relig$ideo3[which(impdat20_relig$politicalid7 < 4)] <- 1
impdat20_relig$ideo3[which(impdat20_relig$politicalid7 == 4)] <- 0
impdat20_relig$ideo3[which(impdat20_relig$politicalid7 > 4)] <- -1
impdat20_relig$ideo3_lab <- factor(impdat20_relig$ideo3,
                                labels = c("liberal", "moderate", "conservative"))

impdat20_relig$ideo_lib <- ifelse(impdat20_relig$ideo3 == -1, 1, 0)
impdat20_relig$ideo_mod <- ifelse(impdat20_relig$ideo3 == 0, 1, 0)
impdat20_relig$ideo_con <- ifelse(impdat20_relig$ideo3 == 1, 1, 0)

impdat20_relig$ideo7 <- impdat20_relig$politicalid7*-1 + 7
impdat20_relig$ideo7_sc <- impdat20_relig$ideo7/6

# Region
impdat20_relig$south <- NA
impdat20_relig$south <- ifelse(impdat20_relig$STATE %in% c("VA", "NC", "SC", "GA", "FL", "MS", "AL", "TN", "AR", "LA", "TX"),
                            1, 0)
impdat20_relig$south[impdat20_relig$STATE == ""] <- NA

impdat20_relig$census_region <- NA
impdat20_relig$census_region[which(impdat20_relig$STATE %in% c("DE", "MD", "DC", "WV", "KY", "OK",
                                                         "VA", "NC", "SC", "GA", "FL", "MS", "AL", 
                                                         "TN", "AR", "LA", "TX"))] <- "South"
impdat20_relig$census_region[which(impdat20_relig$STATE %in% c("CA", "OR", "WA", "MT", "WY",
                                                         "ID", "CO", "UT", "AZ", "NM",
                                                         "NV", "AK", "HI"))] <- "West"
impdat20_relig$census_region[which(impdat20_relig$STATE %in% c("ND", "SD", "MN", "WI",
                                                         "IA", "NE", "KS", "MO",
                                                         "IL", "IN", "MI", "OH"))] <- "Midwest"
impdat20_relig$census_region[which(impdat20_relig$STATE %in% c("ME", "VT", "NH", "MA", "RI",
                                                         "CT", "NY", "NJ", "PA"))] <- "Northeast"
impdat20_relig$census_region[impdat20_relig$STATE == ""] <- NA
impdat20_relig$census_region <- as.factor(impdat20_relig$census_region)


# Education
impdat20_relig$edu_cat <- NA
impdat20_relig$edu_cat[which(impdat20_relig$edu <= 3)] <- 1 # Less than HS
impdat20_relig$edu_cat[which(impdat20_relig$edu == 4)] <- 2 # HS 
impdat20_relig$edu_cat[which(impdat20_relig$edu == 5 | impdat20_relig$edu == 6)] <- 3 # Some college or associates
impdat20_relig$edu_cat[which(impdat20_relig$edu == 7)] <- 4 # College
impdat20_relig$edu_cat[which(impdat20_relig$edu >= 8)] <- 5 # Adv degree or schooling
impdat20_relig$col <- ifelse(impdat20_relig$edu_cat >= 4, 1, 0)

impdat20_relig$edu_cat4 <- NA
impdat20_relig$edu_cat4[which(impdat20_relig$edu_cat <= 2)] <- "HS or Less"
impdat20_relig$edu_cat4[which(impdat20_relig$edu_cat == 3)] <- "Some College"
impdat20_relig$edu_cat4[which(impdat20_relig$edu_cat == 4)] <- "College"
impdat20_relig$edu_cat4[which(impdat20_relig$edu_cat == 5)] <- "Post-Grad"
impdat20_relig$edu_cat4 <- factor(impdat20_relig$edu_cat4,
                               levels = c("HS or Less",
                                          "Some College",
                                          "College",
                                          "Post-Grad"))

# Reason for Visiting
# impdat20_relig$broughtwebsite
impdat20_relig$broughtwebsite2 <- NA
impdat20_relig$broughtwebsite2[which(impdat20_relig$broughtwebsite %in% c("Assignment for school", "Assignment for work"))] <- "Assignment school/work"
impdat20_relig$broughtwebsite2[which(impdat20_relig$broughtwebsite %in% c("Mention in a news story (any medium)"))] <- "News Mention"
impdat20_relig$broughtwebsite2[which(impdat20_relig$broughtwebsite %in% c("My Internet search for this topic or a related topic", "Mention or link at a non-news Internet site"))] <- "Internet Mention/Search"
impdat20_relig$broughtwebsite2[which(impdat20_relig$broughtwebsite %in% "Recommendation of a friend of co-worker")] <- "Peer Recommendation"
impdat20_relig$broughtwebsite2[which(impdat20_relig$broughtwebsite %in% c("Other", "null"))] <- "Other/null"

impdat20_relig$broughtwebsite2 <- as.factor(impdat20_relig$broughtwebsite2)


impdat20_relig$broughtwebsite4 <- NA
impdat20_relig$broughtwebsite4[which(impdat20_relig$broughtwebsite %in% c("Assignment for school", "Assignment for work"))] <- "Assignment school/work"
impdat20_relig$broughtwebsite4[which(impdat20_relig$broughtwebsite %in% c("My Internet search for this topic or a related topic", "Mention or link at a non-news Internet site", "News Mention"))] <- "News/Internet Mention/Search"
impdat20_relig$broughtwebsite4[which(impdat20_relig$broughtwebsite %in% "Recommendation of a friend or co-worker")] <- "Peer Recommendation"
impdat20_relig$broughtwebsite4[which(impdat20_relig$broughtwebsite %in% c("Other", "null"))] <- "Other/null"

impdat20_relig$broughtwebsite4 <- as.factor(impdat20_relig$broughtwebsite4)

# Interview Date Indicators -----------------------------------------------
impdat20_relig$post_gf <- ifelse(impdat20_relig$date_mdy > "2020-05-25", 1, 0)
impdat20_relig$post_gf_1month <- ifelse(impdat20_relig$date_mdy > "2020-05-25" & impdat20_relig$date_mdy <= "2020-06-26", 
                                     1, 0)
impdat20_relig$post_gf_2month <- ifelse(impdat20_relig$date_mdy > "2020-06-26" & impdat20_relig$date_mdy <= "2020-07-26", 
                                     1, 0)
impdat20_relig$post_gf_3month <- ifelse(impdat20_relig$date_mdy > "2020-07-26" & impdat20_relig$date_mdy <= "2020-08-26", 
                                     1, 0)
